Angular dependence of the upper critical field of Sr2Ru04 
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One of the remaining issues concerning thie spin-triplet superconductivity of Sr2Ru04 is the strong limit of 
the in-plane upper critical field //c2 at low temperatures. In this study, we clarified the dependence of H^i on 
the angle 6 between the magnetic field and the ab plane at various temperatures, by precisely and accurately 
controlling the magnetic field direction. We revealed that, although the temperature dependence of H^2 for 
\d\ > 5° is well explained by the orbital pair-breaking effect, Hc2{T) for \d\ < 5° is clearly limited at low 
temperatures. We also revealed that the //(.2 limit for 16*1 < 5° is present not only at low temperatures, but also at 
temperatures close to Tc- These features may provide additional hints for clarifying the origin of the H^2 limit. 
Interestingly, if the anisotropic ratio in Sr2Ru04 is assumed to depend on temperature, the observed angular 
dependence of //c2 is reproduced better at lower temperature with an effective-mass model for an anisotropic 
three-dimensional superconductor. We discuss the observed behavior of H^i based on existing theories. 

PACS numbers: 74.25.Dw, 74.70.Pq 



I. INTRODUCTION 

The layered perovskite superconductor Sr2Ru04 with the 
transition temperature Tc of 1.5 K has been extensively stud- 
ied due to its unconventional pairing state.' - Knight shift 
measurements with NMR^ and with spin-polarized neutron- 
scattering^ have revealed the invariant spin susceptibility 
across for H \\ ab, which firmly indicates that the spin 
part of the Cooper-pair state is triplet. The orbital part is 
favorably interpreted as odd parity based on the measure- 
ments of the critical current 1^ through Pb/Sr2Ru04/Pb prox- 
imity junctions^ and other experiments These results es- 
tablish that Sr2Ru04 is an odd-parity spin-triplet supercon- 
ductor. In addition, the (iSRiSi and Kerr effectii measure- 
ments indicate broken time-reversal symmetry in the super- 
conducting state. The zero-field ground state consistent with 
all these results is expressed by the vector order parameter, the 
c/-vector, d = Aozikx- + iky). However, recent Ru-NMR mea- 
surements under very low fields down to 20 mT revealed the 
invariant Knight shift for H || cJ^ii^ This means d \\ ab {± H) 
with the following two possibilitiesj '^''^'*'* The c/-vector can 
rotate freely in the ab planer^ or the tZ-vector pointing along 
the c axis in zero field can flip perpendicular to the c axis by a 
small magnetic field along the c axis^^ In either case, the spin 
of the Cooper pair can be polarized to any field directions at 
least above 20 mT. 

Another unsolved issue in Sr2Ru04 is the origin of the 
strong limit of the upper critical field //c2, which occurs when 
a magnetic field is applied parallel to the ab plane. - Similar 
Hc2 limit is observed in another spin-triplet superconductor 
UPt3 for H II cr^ as shown in the inset of Fig.|2tb). These //c2 
limits are reminiscent of the Pauli effect, which results from 
the Zeeman energy of quasiparticles. However, in spin-triplet 
superconductors, the Pauli effect contributes to pair-breaking 
only when d \\ H, because the spin of the triplet Cooper pairs 
can be polarized along the field direction when d ± H. As 
mentioned above, the c/-vector of Sr2Ru04 is likely to be per- 
pendicular to the magnetic field possibly except at low fields. 
This suggests that the Pauli effect should not affect H^t- We 
note that the c/- vector of UPt3 for H \\ c was revealed to be per- 



pendicular to the magnetic field (d || a) in phase There- 
fore, the Hc2 limit observed in UPt3 cannot be attributed to 
the Pauli effect either. The origins of these Hc2 limits have not 
been clarified yet. 

In this paper, we report the dependence of H^i of Sr2Ru04, 
determined from the ac susceptibility, on the magnetic field di- 
rection between the ab plane and c axis. Because of the large 
anisotropy of //c2 in Sr2Ru04, a small misalignment would 
lead to a large difference in the value of //c2, especially when 
the field direction is nearly parallel to the ab plane. Therefore, 
in this study, we controlled the applied field direction more 
accurately and precisely than in the previous reportsiaiii^ We 
evaluated the H^iiT) curves for different field directions and 
revealed that the H^i limit is clearly observed only when the 
angle between the magnetic field and the ab plane is less 
than 5 degrees. This limit was revealed to occur not only 
at low temperatures, but also at temperatures close to T^- We 
also identified the angle 6 dependence of Hc2 at several fixed 
temperatures. We found that Hc2{0) is fitted better at lower 
temperature with an effective-mass model for an anisotropic 
three-dimensional superconductor In addition, we investi- 
gated the difference between Hc2{0) for fields in the (100) 
plane and for fields in the (110) plane. The difference ap- 
pears only at low temperatures below roughly 1 K for small 0. 
These results allow us to reexamine the origin of the Hc2 limit 
based on existing theories. 



n. EXPERIMENTAL 

We used single crystals of Sr2Ru04 grown by a floating 
zone methodisi In this paper, we focus on the result obtained 
from a single crystal with dimensions of approximately 1 .0 x 
0.5 mm- in the ab plane and 0.08 mm along the c axis. The 
directions of the tetragonal crystallographic axes of the sam- 
ple were determined from x-ray Laue pictures. We shaped 
the sample so that the side surface of the sample was 10 
degrees away from the (100) plane in order to avoid possi- 
ble anisotropy effects due to surface superconductivity f22 The 
crystal was annealed in oxygen at 1 atm and 1050 °C for a 
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FIG. 1: (Color online) Field dependence of Ax' = ^'(0.9 K)-y(2 K) 
at several 6 for (a) = 0° and (b) if> = 45°. Here, if) denotes the 
azimuthal angle within the ab plane between the magnetic field and 
the [100] axis and 6 denotes the angle between the magnetic field 
and the ab plane. H^i is defined as the intersection of the linear 
extrapolations in Ax'- Thick and thin arrows represent H^i and the 
anomaly due to the mosaic structure, respectively. The dip in Ax' 
near //;.2 is attributable to the ordinary peak effect— 



week to reduce the amount of oxygen deficiencies and lat- 
tice defects. A sharp superconducting transition was observed 
in the ac susceptibility measurements with the midpoint at 
rc=1.503 K. 

We measured the ac magnetic susceptibility Xac - x' ~ (y" 
by a mutual-inductance technique using a lock-in amplifier 
with a frequency of 887 Hz. The sample was cooled down 
to 70 mK with a ''He-'*He dilution refrigerator The ac mag- 
netic field of 20 ^iT-rms was applied nearly parallel to the c 
axis with a small coil. The dc magnetic field was applied us- 
ing the "Vector Magnet" system^l with which we can con- 
trol the field direction three dimensionally and precisely. The 
accuracy and precision of the field alignment with respect to 
the ah plane are better than 0.1 degree and 0.01 degree, re- 
spectively. Owing to the high sensitivity of the pick-up-coil, 
parasitic background contributes to the x-^c signal.'^A In order 
to obtain ;^^ac contribution only from a superconductivity A^', 
we adopt Ay(r,//) = y(r,//) -^(2 K,//). The small de- 
viation of the normal state values from zero indicates a good 
reliability of the background subtraction. We define H^o as 
the intersection between the linear extrapolations of Ay'(//) in 
the superconducting and normal states, as illustrated in Fig.[T] 
with dashed lines. The directions of the crystalline axes [100] 
and [001] with respect to the field direction were calibrated 
by making use of the anisotropy in //(;2.~^'~' Our highly accu- 
rate and precise measurements revealed that the present sam- 
ple has a mosaic structure dominated by two domains sharing 
the [100] axis; the [001] axis of one domain is tilted nearly 
toward the [010] axis by 0.5 degree from the [001] axis of the 
other part. 
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FIG. 2: (Color online) (a) Field-temperature {H - T) phase diagram 
of Sr2Ru04 at = 0°, 0.5°, 1°, 2°, 5°, 10°, and 90° from top to bottom 
for = 0°. (b) Temperature dependence of h' defined as eq. (Ujl. The 
inset represents the H-T phase diagram of UPtj for// II c(Ref.[ll). 
The dashed curve in (b) is h*{t) of the boundary of the ( B + C ) phase 
in UPt3 for H \\ c (the dashed curve in the inset). 



III. RESULTS 

Figure |2ja) is the field-temperature {H - T) phase diagram 
in various field directions at (p - 0°, where (p denotes the az- 
imuthal angle within the ab plane between the magnetic field 
and the [100] axis. At = 0°, no anomaly due to the mosaic 
structure was seen in the raw Xac data, as shown in Fig. [TJa)- 
Reflecting the large anisotropy of H^i in Sr2Ru04, Hc2 be- 
comes rapidly small when the angle between the magnetic 
field and the ab plane increases from 0°. In the specific heat 
measurements, the second superconducting transition was ob- 
served just below //c2 at low temperatures below 0.8 KM Al- 
though such an additional transition was observed below 0.6 
K in the ac susceptibility measurements, it was difficult to un- 
ambiguously identify it to be attributable to the second super- 
conducting transition.?' This is also the case for the present 
study. One possible reason for this difficulty is that ac sus- 
ceptibility, mainly probing the vortex movements, may not be 
sensitive to the small change in the entropy detected by spe- 
cific heat measurements. Therefore, we do not focus on the 



3 



< 

a: 

0) -0.5 

_N 

IB 
E 

o 

c 

I- -1 

< 

< 



, UPt3(H//c) 
□\N, \ ~-x\x \ 

90° "v,''*-vo:5,k''"^\ 



Sr2Ru04 



0.1 K 



0.5 

t= T/T^ 




30 

d (deg) 



FIG. 3: (Color online) Temperature dependence of the slope of the 
H-T phase diagram of Sr2Ru04 at = 90°, 5°, 2°, r,0.5°, and 0° 
from bottom to top for = 0°. The slope of the //c2 curve for UPts 
(Refill), evaluated from the dashed curve in the inset of Fig.[2jb), 
is also plotted with crosses. The dashed curves are guides to the eye. 



feature of the additional transition in this paper 

To characterize the limit of H^iiT), we normalized H^i by 
the initial slope at T^: 



h\f) = - 



d//c2/df|,=i 



(t = TIT,). 



(1) 



If Hc2 is determined by the orbital pair-breaking effect, 
which originates from the kinetic energy of supercurrent 
around magnetic vortices, i/c2 is described by the Werthamer- 
Helfand-Hohenberg (WHH) theor y ^^-^^ and its extension to p- 
wave superconductors In these theories, it is expected that 
h*(f) increases linearly on cooling and is weakly suppressed at 
low temperatures with h*{t = 0) ~ 0.7. In Fig.|2b), we plot 
h*{f) with different 9at(p - 0°. The initial slope d//c2/df|,=i is 
defined from the linear fit to Hdit) in the region 0.85 < f < 1. 
For l^l > 5°, li*(t) behaves as expected from the WHH theory. 
In contrast, for \6\ < 2°, h*(t) is strongly limited at low temper- 
atures. This result indicates that the Hd limit in Sr2Ru04 is 
prominent for \d\ < 5°. To emphasize the //c2 limit in another 
spin-triplet superconductor UPt3, we plot, in Fig.|2jb) with the 
dashed curve, h*{t) of the boundary of the ( B -H C ) phase for 
H \\ c. Although we chose the less limited one between the 
two Hc2{T) curves in UPt3 for H \\ c, a strong Umit of li*(t) is 
clearly seen. 

If the orbital pair-breaking effect is mainly responsible for 
determining Hd, the slope of the H-T phase diagram should 
be constant down to well below Tc. To identify the limit of 
Hc2 near T,, we evaluate the slope at temperature {Ti + T2)/2 
as AHd/AT = [H^iTi) - H,2(T2)]/(Ti - T2), where Ti and 
T2 are temperatures of adjacent data points. The results are 
shown in Fig. [3] For \0\ > 5°, the slope is constant down to 
approximately 1 K and approaches zero at low temperatures, 
which is well explained by the orbital pair-breaking effect. 
However, for \6\ < 5°, the slope near is not temperature- 
independent any more. This result suggests that the Umit 
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FIG. 4: (Color online) (a) Field-angle 6 dependence of the upper 
critical field H^i of Sr2Ru04 at various temperatures for ip = 0° . (b) 
Angle 9 dependence of H^2 normalized by Hc2(S = 90°). Solid curves 
are the fitting results using eq. ([3} in the range 2.5° < \d\ < 90°. Insets 
are enlarged views near 9 = 0°. 



observed for \0\ < 5° is present not only at low temperatures, 
but also at temperatures close to T,. The slope of the H^iT) 
curve of UPt3 for H \\ c (the dashed curve in Fig. Hfb)) also 
continues to vary up to T,, as plotted in Fig. [3] 

In order to characterize the non-linear temperature depen- 
dence of Hc2 in Sr2Ru04, we fitted //c2(0 for 6 - 0° hy 
a{l - t)" with fitting parameters a and n. In any fitting range, 
n is obviously larger than n - 0.5, which is expected for the 
two-dimensional (2D) superconductivity;.^ the fitting in the 
range 0.9 < f < 1 yields n = 0.9. In addition, the coherence 
length along the c axis is estimated to be 3.2 nm using the 
GL equation 



(2) 



with HQHc2\\ab = 1.5 T and jioHdWc - 0.075 T. Here, <t>{) is 
the flux quantum. This value of is five times larger than 
the spacing of the conductive RUO2 layers (0.62 nm)?^ Even 
if the WHH value -0.7d//c2/df|,=i = 2.5 T is used for H^Wab, 
= 1.9 nm is obtained. These facts indicate that the super- 
conductivity of Sr2Ru04 cannot be classified as a 2D super- 
conductivity. 

Figure |4ja) represents the dependence of Hd at various 
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FIG. 5: (Color online) Error of the fitting of eq. ([3} to the H^2{0) data 
for = 0° in the range 6(n,i„ < 16*1 < 90° for several d,„,„. 0niin = 2.5° 
is likely to give the most appropriate fitting range. 



temperatures for = 0°. The 6 dependence of //c2 normalized 
by Hc2(0 = 90°) is also plotted in Fig. Etb). We found that, 
although H^2{e)/H,2{90°) for 10° < \e\ < 90° is nearly inde- 
pendent of temperature, it decreases on cooling for \6\ < 10°. 

We found that the observed Hc2{ff) is well explained 
by the Ginzburg-Landau (GL) theory for anisotropic three- 
dimensional (3D) superconductors,— if we allow the 
anisotropic ratio F - //c2(0°)///c2(90°) to depend on temper- 
ature. The angular dependence of //c2 is expressed as 



^nry- sin- e + cos^ e 



(3) 



We fit eq. ^ to the observed Hc2{ff) at temperature Tq with 

jorb ( 



two fitting parameters H°^f^^i^(TQ) and r(To). We chose the fit- 



TABLE I: The upper critical field and its anisotropy. Fitting param- 
eters are obtained by the fit of H^2iff) for ^ = 0° using eq. l[3j with 
'^min = 2.5° at each temperature. 
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FIG. 6: (Color online) Comparison between Ha(d) al (p = 0° (cir- 
cles) and = 45° (crosses). 



ting range as 9^i„ < \6\ < 90° so that the range is as wide 
as possible while the fitting yields a good result in the whole 
chosen range. Figure |5] represents the error of the fitting for 
different d^i„ at each temperature. When O^in > 2.5°, Hc2{8) is 
well fitted by eq. (O in the chosen fitting range. By contrast, 
for flniin < 2.5°, Hc2{0) exhibits systematic deviation from 
eq. ([3]) around \0\ ~ 2°. Thus, we conclude that 0^i„ = 2.5° is 
the most appropriate. The fitting results are plotted in Fig. |4] 
with the solid curves and the obtained fitting parameters are 
listed in Table U Interestingly, the observed Hc2{d) is fitted by 
eq. (O better at lower temperatures, as being clear in the in- 
set of Fig.HJb). This tendency is also clear when the fit ratio 
r is compared with the experimental ratio //c2(0°)///c2(90°). 
We should mention that the thin-film model^^ applied to a 2D 
superconductor^^ in which Hc2{d) exhibits a cusp at = 0°, 
cannot account for our data. While we carefully examined the 
9 dependence of Hc2, a kink in Hc2{0) around 0-2° revealed 
by the specific heat measurements at 0.1 K (Ref . [l7h was not 
detected in the present study. We note that a kink in Hc2{0) 
was not detected in the thermal conductivity measurements 
at 0.32 K, either (Fig. 4(a) in Ref. [H). On the basis of the 
presently available results, we cannot clarify why the kink in 
Hc2{0) was observed only in the specific heat measurement at 
0.1 K. 

In Fig. |6] we compare the dependence of i/c2 at angles 
between = 0° and (p - 45°. As indicated by Fig. [lib). 
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FIG. 7: (Color online) Temperature dependence of the in-plane //c2 
anisotropy between = 0° and ip = 45° at = 0°. A(ji(,H^2) is 
defined as noH^ii^ = 45°) - /jo^c2(0 = 0°). 



two onset features appear in the field dependence of Ax' at 
(p - 45°, reflecting the fact that the present sample consists 
mainly of two domains. Since it is possible to separate the 
contribution from each of these domains, we plot in Fig |6] 
Hciiff) for the major domain. From Fig. |6] we found that 
HciiT, 6) at <p - 45° is both qualitatively and quantitatively 
similar to Hc2{T,6) at 4> - 0°. Smafl diff'erence in H^iiO) at 
angles between = 0° and = 45° is observed only when the 
magnetic field is applied nearly parallel to the ab plane. Fig- 
ure |2] represents the temperature dependence of the in-plane 
Hc2 anisotropy between = 0° and (p - 45°. Here, we define 
A(azo//c2) as ^ioH,2{4> = 45°) - ^iqH,2((P = 0°). Although the 
temperature T* at which A{fioHc2) starts to increase on cool- 
ing depends on samples (0.9 K < T* < 1.2 K), A(/io//c2) of 
40 mT at low temperatures is nearly the same among diff'erent 
samples with best T^i^ 

On the basis of the phenomenological theory proposed by 
GorkoVj-4 superconductivity with a two-component order pa- 
rameter, k f + iky, should be accompanied by a substantial four- 
fold anisotropy in the in-plane Hc2- However, as presented in 
Fig. 121 no in-plane Hc2 anisotropy is observable above about 
1 K; the anisotropy grows on cooling, but reaches at most 
3% at low temperatures. This lack of the large in-plane //c2 
anisotropy is attributable to the multiband effect] -^^-^^-^^ Be- 
cause the directions of the gap minima are 45 degrees differ- 
ent between the active (y) and passive (a and /3) bands^^ the 
Hc2 anisotropy reflecting the gap structure on different Fermi 
surface sheets can be cancelledi^^i^ 



Sr2Ru04 is an unusual orbital pair-breaking effect. For ex- 
ample, in a nearly 2D superconductor (TMET-STF)2BF4i^ it 
is proposed that Hc2{T) for H \\ a is limited due to the limit 
of the coherence length by the layer spacing, which leads to 
the decrease of the anisotropy ratio of Hc2 on cooling. For 
Sr2Ru04, estimated using eq. for which the ordinary or- 
bital pair-breaking effect is assumed, is limited to be approxi- 
mately 3.2 nm below 1 K. In fact, if we strictly apply eq. (|2|, 

takes a minimum at 0.8 K and even increases by about 3% 
at low temperatures. However, we cannot find a clear answer 
to this limit because the limited value of 3.2 nm, is five 
times larger than the layer spacing. Therefore, the origin of 
the apparent limit of the coherence length in Sr2Ru04 seems 
different from that in (TMET-STF)2BF4. 

Recently, Machida and Ichioka proposed the Pauli effect as 
an additional pair-breaking effect leading to the Hc2 limit in 
Sr2Ru04.— Using a model with a single-band spherical Fermi 
surface and by assuming the Pauli effect, they reproduced the 
observed H^2{d) at 0. 1 K (Ref.[il) as well as field dependences 
of the specific heat^ and magnetization.— Interestingly, we 
found that the Machida-Ichioka model well reproduces our 
results of Hc2{T, 0), too. Nevertheless, this would not lead to 
the conclusion that the Hc2 limit in Sr2Ru04 is attributable 
to the Pauli effect because the Machida-Ichioka model over- 
looks some key experimental as well as theoretical facts. First, 
Machida-Ichioka model does not include the multiband ef- 
fect. Their single-band model explains the field dependence 
of the specific heat at low temperatures. However, Sr2Ru04 
has three cylindrical Fermi surfaces, a, f5, and yi^^ Although 
the active band y is dominant in the superconductivity in high 
fields, the passive bands a and /3 also contribute to the super- 
conductivity in low fields.— The contribution from the a and 
/3 bands is essential to explain the plateau-like dependence 
quantitatively^^ In fact, inclusion of the multiband effect is 
needed to explain the dependence of the specific heat at 
low temperatures in zero field-^^^*^^^ Secondly, as mentioned 
in Sec. I, the Pauli effect contradicts the results of the Knight 
shift experiments jMjS, Although they proposed the possibility 
that the spin part of the Knight shift was too small to be de- 
tected in the NMR experiments, the spin part at the Ru site is 
in reality as large as 4%,'* In addition, the superconductivity 
was distinctly observed in l/Ti through Ru NMR in the iden- 
tical setup.— These facts exclude the possibility of the Pauli 
mechanism. Therefore, an alternative mechanism needs to be 
introduced to explain both the Knight shift behavior and the 
Hc2 limit. 



V. SUMMARY 



IV. DISCUSSION 

Let us discuss the origin of the Hc2 limit in Sr2Ru04. For 
both 2D and 3D superconductors in which the main pair- 
breaking effect is due to the ordinary orbital effect, such a 
limit is not expected. Thus, in order to explain the Hc2 limit, 
we need an additional pair-breaking mechanism. 

One of the possible additional pair-breaking effects in 



We have clarified the temperature and field-angle depen- 
dence of Hc2 of Sr2Ru04. Our experiments were performed 
with an accurate and precise control of the applied magnetic 
field to avoid errors due to the misalignment. We revealed 
that the Hc2 limit is clearly observed for \9\ < 5° and it occurs 
not only at low temperatures but also at temperatures close to 
Tc- We also found that, by assuming a temperature-dependent 
anisotropic ratio, the GL theory for an anisotropic 3D super- 



6 



conductor can explain the angular dependence of H^2 well, 
particularly at lower temperatures. The observed behavior of 
Hciiff) is qualitatively the same between <p - 0° and cp = 45°. 
Only a small in-plane H^i anisotropy was observed at low 
temperatures, which disappears rapidly as the magnetic-field 
direction leaves from the ab plane. Until now, the origin of the 
effective pair-breaking effect, which is compatible with both 
the invariance to the Knight shift and the hmiting behavior of 
Hc2, remains unclear. 
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